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Abstract 

In this paper we consider the collision of black holes with parallel spins using first order pertur- 
bation theory of rotating black holes (Teukolsky formalism). The black holes are assumed to be 
close to each other, initially non boosted and spinning slowly. We estimate the properties of the 
gravitational radiation released from such an collision. The same problem was studied recently by 
Gleiser et al. in the context of the Zerilli perturbation formalism and our results for waveforms, 
energy and angular momentum radiated agree very well with the results presented in that work. 
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I. INTRODUCTION 



There is considerable current interest in studying the collision of two black holes, since 
these events could be primary sources of gravitational waves for interferometric gravitational 
wave detectors currently under construction. The close limit approximation applies to the 
late stage of such an event, when the system can be considered as a single distorted black 
hole. One can then use linear perturbation theory to estimate the energy and angular 
momentum lost by the system due to gravity wave emission, and even obtain "waveforms" 
which could be useful for experimental detection of these waves. 

In the past, this method has been applied to the case of a head-on collision of two boosted 
holes and the case of slow inspiral @ with considerable success. Recently, this method was 
extended to cover the collision of spinning holes with anti-parallel spin || and case of equal 
and parallel spin [Q. However in the latter work, the close limit of the two merging holes 
was considered as a distorted Schwartzchild hole. Therefore, the perturbation formalism 
used in that work, was the Zerilli formalism. In this paper, we shall treat the merger of two 
equal mass holes with equal and parallel spin, as a distorted Kerr black hole. Thus, we shall 
use the Teukolsky formalism for our evolutions, and compare our results with those from 
the Zerilli formalism. 

We also hope that such a comparative study, will shed light on the apparent discrepancy 
that was noted for the amount of the angular momentum lost for the case of slow inspiral 
of two equal mass holes, as obtained by these two perturbation formalisms 0. Comments 
relating to that shall be published elsewhere. 

It should be noted that, results presented here can be easily combined with results ob- 
tained in the past, (say) for the slow inspiral case H via simple superposition, to obtain 
waveforms, etc. for an event in which two equal mass and equal and parallel spin holes 
merge. This is because all these results are based on linear black hole perturbation theory . 

II. INITIAL DATA 

To evolve a spacetime in general relativity, one needs as initial data, a 3-geometry g a t, and 
an extrinsic curvature K ab , that solve Einstein's equations on some starting hypersurface. 
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These initial value equations have the form, 

V a (K ab -g ab K) = (1) 
3 R — K ab K ab + K 2 = (2) 

where 3 R is the scalar curvature of the three metric. If we propose a 3-metric that is 
conformally flat g ab = <p A g a b, with g ab a flat metric, and </> 4 the conformal factor, and we 
use a decomposition of the extrinsic curvature K ab = (fi~ 2 K ab , and assume maximal slicing 
K% = 0, the constraints become [§], 

V a K ab = (3) 
V 2 = -h- 7 K ab K ab , (4) 

o 

where V is a flat-space covariant derivative. 

To solve the momentum constraint, we start with the well known solution that represents 
a single hole with spin S, 
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e C adS n c n b + e cbd S d n c n a . (5) 



In this expression for the conformally related extrinsic curvature at some point x a , the quan- 
tity n b is a unit vector, in the flat space with metric g ab , directed from a point representing 
the location of the hole to the point x a . The symbol R represents the distance, in the flat 
base space, from the point of the hole to x a . It is straightforward to show that the solution of 
the Hamiltonian constraint corresponding to equation (5) corresponds to a spacetime with 
ADM angular momentum S a . 

The next step is to modify this to represent holes centered at x = ±L/2 in the flat metric. 
Since the momentum constraint is linear, we can simply add two expressions of the above 
form, 

KT = KT (x ^ x ~ L/2,S Z = S) + KT (x ^ x + L/2,S Z = S) . (6) 

We will now use a polar coordinate system in the flat space determined by g ab centered in 
the mid-point separating the two holes and label the polar coordinates as (R,9,(p). Thus, 
R will be the distance in the flat space from the midpoint between the holes. 

To solve the Hamiltonian constraint |], we use the "punctures" anzatz, i.e. we assume 
that is of the form, 

<t> = 4>BL + 4>Reg (7) 



where, 

M M 

^ = 1 + 2^ + 2^ (8) 
that is, 4>bl is taken as the Brill-Lindquist conformal factor ||, and we demand that (j)R eg 
must be regular in the whole conformal plane, and vanish for large R{. Here, R\ and R2 
indicate the distance measured in conformal space, from the each of the two holes. If <ft is 
to be a solution of equation |], (pR eg must satisfy the equation, 

1 K ab K i. 



where K ab is equal to K^°, as defined in equation ^. 

It is straightforward to see that the total initial K ab scales linearly with S. If we assume 
that S and L are small quantites (the holes are close to each other initially, and spinning 
slowly) and are of the same order, we need to keep only the quadratic terms in S in the 
the source term of the Hamiltonian constraint, and we can neglect the rest. Then the 
Hamiltonian constraint would become, 

v ^ = -»WiF (10) 

where, J = 2S the total angular momentum of the system. We solve this equation, following 
Gleiser et al. [|j] and obtain the final result for accurate to order J 2 , 

<Pne 9 = <j> m {R)Y2(9,<p) + <j> m {R)Y°(9,<p) (11) 



where 



40FJ 2 (M 4 + 10RM 3 + A0R 2 M 2 + 40R 3 M + 16i? 4 ) 
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5M 3 (2R + M) 

, D , 32 v / o7j 2 i? 2 , in , 

$(2 ' 0)(i?) = -2hM{2R + Mf (12) 

For the sake of completion we also list here explicitly, the components of extrinsic curvature, 
keeping only the lowest two orders, 

. JX 2 sin(2^) sin 2 (#) 



Krr = - 
Kr B = - 

K Rip = 3 



R 5 

9 JL 2 sin(2y?) sin(0) cos(#) 
"8 ~R A 

Jsin 2 (g) _ 3 JL 2 (7 + cos 2 (y?) (25 cos 2 (fl) - 19))sin 2 (fl) 
~R 2 8 W 



K ee 



-JL 2 sin(2^) cos 2 (6)/R 3 



K, 




(13) 



We must now map the coordinates of the initial value solution to the coordinates for 
a Kerr black hole with a = J/M background. To do this, we interpret the R as the 
isotropic radial coordinate and we relate it to the usual Boyer-Linquist radial coordinate r 
by R = (y/r + s/r - 2M) 2 /4. From this we arrive at the final perturbative expressions for 
the components of the metric and extrinsic curvature. We shall interpret these quantities 
as perturbations over a Kerr solution background with perturbation parameter L 2 . It is 
important to note that we computed the initial data only to order J 2 mainly for ease in 
analytic computation, and we do not interpret J as a formal perturbation parameter. Since 
we are using the Teukolsky formalism for our evolutions, J is simply a background quantity 
and not a perturbation parameter. This is an important difference from the corresponding 
Zerilli computation, where J is a perturbation parameter on equal footing with L. 

Using those expressions for the metric and extrinsic curvature, we calculate initial data for 
the Teukolsky function \1/ following the prescription provided in reference . The expressions 
we arrived at, are far too long and complicated to include in this paper, and to be of any 
direct use. Explicit algebraic expressions for this initial data shall be provided (in the 
computer algebra system, MAPLE format or as FORTRAN code), on request, to anyone 
interested. 

III. EVOLUTION OF THE DATA USING THE TEUKOLSKY EQUATION 

Given the Cauchy data from the last section, the time evolution is obtained from the 
Teukolsky equation ||, 




M(r 2 - a 2 ) 
A 



(15) 



(14) 
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where M is the mass of the black hole, a its angular momentum per unit mass, E = 
r 2 + a 2 cos 2 9, and A = r 2 — 2Mr + a 2 . Evolving the initial data we just calculated, with this 
equation will enable us to extract gravity wave waveforms that correspond to the late stage 
merger of spinning holes. We use the 2 + 1 dimensional Teukolsky evolution code written by 
Krivan et al. || using a radial (tortoise coordinate grid) resolution of M/20 and an angular 
resolution of ir/AO. We can also estimate the energy carried away by these gravitational 
waves using |10| . 



d -^= lim/JL / dn 



dt r ^°° 4nr 6 Jn 



dt V(t,r,9,(p) 



dQ = sin 9 d-d dip. (16) 



The angular momentum radiated can similarly be calculated using [10 

ft \ ( rt rt' 



= — lim | —^—z Re 

dt r ^°° [A7rr b 



f dn(d v [ diV(i,r,e,<p))( f dt' [ diV(i,r,6,<p) 

JfTl \ J— oo / \J— oo J —oo 



(17) 



IV. RESULTS OF THE EVOLUTIONS 



In this section we show waveforms and plots for energy and angular momentum radiated 
from the collision of two spinning holes (with parallel spin) and compare the results to those 
obtained from the Zerilli formalism. Recall that the two holes have equal mass and equal 
and parallel spin. 

The waveforms that follow, are for a collision of two black holes that were initially sep- 
arated by a conformal distance of 1.0 and have an individual spin of 0.05 in units of ADM 
mass. The waves were extracted at radial location with r* = 25M and at a polar angle 
9 = tt/2. 

In figure 1 we show the m = 2 mode of the Teukolsky function as a function of time. We 
see the typical quasi-normal ringing, in both the real and imaginary parts of the function. 
Note that the imaginary part the of waveform, there appears to be a mixing of frequencies. 
This is exactly what was observed by Gleiser et al. ||. They noted a mixing of £ = 2 and 
i = 3 spherical harmonic modes in their evolutions. This is the type of signal that gravity 



wave observatories like LIGO pi] , will detect if they happen to witness a collision of the 
kind we are considering. 

Let us now turn to the results for the radiated energies. Figure 2 shows the radiated 
energy as a function of the initial total spin, for a fixed separation of the holes. Note that the 



Zerilli results agree very well upto about J/M 2 = 0.2. Beyond that, both these calculations 
cannot expect to yield accurate results. This is because both the methods are based on the 
close, slow approximation. Similarly, our results for radiated angular momentum, figure 3, 
agree very well to about J/M 2 = 0.2 but then diverge beyond that. 

One may make the observation that in both these cases (for large a) the Teukolsky 
formalism results appear to have more radiation than the corresponding Zerilli based results. 
This can be explained, in part, by the fact that there is less damping in the QN modes of a 
Kerr hole. However, it is not clear that this physical reasoning accounts for the larger values 
completely. One must keep in mind, that for larger values of a, both these perturbative 
calculations break down, and therefore it is difficult to make any statements apart from 
observing some rather general trends. 

Also worth noting is that our results suggest that such a collision is unlikely to change 
the inspiral based estimate of 1% of the total system mass being radiated away by 
gravitational radiation. This is because, for every value of J (especially the larger values), 
this collision radiates much less than the inspiral case. One would therefore conclude that 
the radiation from the collision of two black holes is dominated by radiation coming from 
the "inspiral" part 0. This fact was also observed in the context of holes with anti-parallel 
spins ||. There also appears to be no appreciable change in the estimate for the radiated 
angular momentum. 

V. A "REALISTIC" MODEL 

In this work we have treated the collision of two black holes, with like masses and spins. A 
more astrophysically likely event would be one that also has some orbital angular momentum, 
i.e. in addition to the two holes spinning, they are also inspiralling into each other. In this 
section we shall attempt such an evolution. 

As mentioned before in this work, our approach here is going to be to use simple su- 
perposition of waveforms from our earlier work |2| with the ones presented here. We can 
do this, since all these results are based on first-order perturbation theory of black holes. 
In addition to obtaining waveforms, we shall also obtain amounts for energy and angular 
momentum radiated. We choose two equal mass black holes located at x — ±0.4, with equal 
and parallel spins of 0.1 each. The spins are aligned along the z-axis. Next we boost the 
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black holes along the positive and negative y-axis respectively, so as to obtain a total orbital 
angular momentum of 0.2. Thus, the total angular momentum amounts to 0.4. Please note 
that all the above mentioned quantities are in units of ADM mass. 

In Figure 4 we include waveforms from such an evolution. The figure depicts the real 
and imaginary parts of the m = 2 mode of the Teukolsky function. It should be noted that 
this waveform is very close to the waveform produced by the "inspiral" part of this collision, 
indicating that, that is dominant part of such an event. We also obtain estimates for energy 
and angular momentum carried away by these waves, 

E RAD /M ADM = 1.0xlO- 3 

Jrad/M adm 2 = 2.2x10- 3 . 

Moreover, by taking a difference between the energies obtained for a case in which the spins 
are aligned with the orbital angular momentum and the case in which they are anti-aligned, 
we can even estimate the lowest order coupling term between the orbital part and spin part. 
Also, by examining the analytic expressions for initial data, one can easily see that the form 
of this term has to be of the kind, JSL 2 . Here J is the orbital angular momentum, S is the 
total spin and L is the separation between the two holes. Numerical evolutions confirm the 
above and yield an empirical estimate for this term, 

Ecplg/M adm = (1.6xlO- 3 )JSL 2 /M ADM 6 . 

VI. CONCLUSIONS 

We performed a first order, perturbative calculation to study the merger of two spinning 
holes with parallel spin based on the Teukolsky formalism. Our results agree very well with 
those obtained in the same context using the Zerilli formalism. 

We also noted that our results here do not significantly affect the original inspiral based 
estimates for energy and angular momentum radiated. 
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FIG. 1: Real and imaginary parts of the Teukolsky waveform extracted at r* = 25M. All 
quantities are in units of ADM mass. 
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FIG. 2: Radiated energies as a function of total initial angular momentum extracted at r* = 25M. 
Note the good agreement between the two formalisms upto J/M 2 = 0.2. All quantities are in units 
of ADM mass. 
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FIG. 3: Radiated angular momenta as a function of total initial angular momentum r* = 25M. 
Note the good agreement between the two formalisms upto J/M 2 = 0.2. All quantities are in units 
of ADM mass. 
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FIG. 4: Real and imaginary parts of the m = 2 mode of the Teukolsky function. All quantities are 
in units of ADM mass. 
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